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THERMAL EFFECTS ON THE CLEARANCE AND STIFFNESS 
OF FOIL JOURNAL BEARINGS FOR A BRAYTON CYCLE 
TURBOALTERNATOR 


SUMMARY 


An analysis of foil journal bearings for a NASA Brayton Cycle 
Unit (BRU) is presented. The study represents an extension of pre- 
vious work in that it includes the effects of thermal expansion of foil- 
bearing components, as well as an improved model of the influence of 
foil flexure. 

The BRU foil-bearing design, based on parameters previously 
established through a simplified thermal calculation, has been tested 
by the computer program developed in this work. The results presented 
herein estimate the range of the bearing film thickness, the bearing 
stiffness and the foil tension as functions of the expected operating 
temperatures and the elastohydrodynamic and geometrical parameters 
selected in this design. The parameter values chosen in the design to 
meet the requirements of stiffness and gap are confirmed. Results for 
alternative foil material and thickness combinations are also given to 
illustrate off-design conditions. 

It is shown that the foil material is a sensitive variable, which 
the designer has at his disposal for controlling the bearing characteristics 
when other parameters nave Deen firmly established. 
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1.0 INTRODUCTION 


The three principal attributes of foil bearings — (a) stability, 

(b) ability to accommodate thermal distortions, and (c) tolerance of foreign 
particles in the gas film — have been demonstrated in the course of pre- 
vious investigations [l] to [ 7] . The development has progressed to the 
point, at which pivoted-shoe bearings in one experimental Brayton Cycle 
Unit (Fig. 1) are to be replaced with foil bearings of the present design. 

The objective of the analysis and data presented herewith is to 
provide a design tool and means of estimating the performance of simi- 
lar types of foil bearings in smaller, or larger turbomachines. 

Because of the complexities of foil bearing analysis, it has not 
been feasible to combine the time-dependence of the fluid film with a full 
treatment of the thermal, elastic and geometric aspects of the problem. In 
reference L 5] , however, it is shown for the case of a similar but idealized 
configuration, that damping assumes no negative values and that the foil 
bearing is, therefore, stable. At the same time, the stiffness does not 
differ appreciably from the quasi-static value in a wide band of frequencies 
of excitation. Therefore, without degrading the results, the more tractable 
quasi-static approach permits the inclusion of realistic thermal, elastic 
and geometric aspects of foil bearings in the analysis. 

The present work represents an addition and a refinement of anal- 
yses included in references L 1 ] and [6] , in that a more realistic model is 
used in accounting for thermal effects , and refinements are introduced in 
estimating the influence of foil flexure. For continuity of presentation, 
certain portions of these analyses are included here with the added improve- 
ments . 
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PIVOTED -SHOE h. 
BEARINGS TO BE / 
REPLACED BY FOIL L 
BEARINGS \ 



NASA SRAVTON ROTATING UNIT 

[ es ;~| msunaiMMjfAcnimc 

DIVISION OF ARIZONA 


NASA Brayton Cycle Unit with Pivoted Bearings. 
(Manufactured by AiResearch Co.) 

The pivoted-shoe bearings in one unit are being 
replaced by foil bearings designed and manu- 
factured by the Ampex Corporation. 


2.0 FORMULATION 


Consider three "infinitely wide" foil sectors spaced at equal 
angular intervals and wrapped around a journal as shown schematically 
in Fig. 2. Each foil sector is supported by a kidney-shaped "shield- 
damper" and wrapped around the cylindrically milled ends, referred to 
as "guides". The shield-dampers are rigidly connected to plates 
parallel to the plane of the paper. 


In the absence of radial load and fluid film and with a uniform 
temperature / t 0 and preload tension per unit width T 0 in each foil sector, 
the journal center is located at the point O, which will be ta v en as the 
origin. This state will be termed the "reference state". When rotation 
ensues, fluid films build up between the foil sectors and the journal, but 
under steady-state conditions, thermal axial symmetry and no radial load, 
the rotor center remains at O. Although the results of this report are 
limited to these conditions, the evaluation of stiffness requires consider- 
ation of rotor displacements from O. Consequently, the formulation starts 
with a general rotor position and is later specialized for axial symmetry. 

Assume that at some instant, the journal is located at a point (x,y) 
in terms of a fixed Cartesian system of coordinates with origin at O. Con- 
sider, further, a set of three auxiliary coordinate systems x k , y k (k = 
1,2,3) with origins at O, such that y k is directed along the bisector of 

the kth foil arc, positive away from the foil and making an angle 

with the y axis. For the given shaft displacement (x,y) , the coordinates 

(x k ,y k ) are given by the transformation 


Xj,. = K Cos ^ -v vj sin 

+ Jj Cos 
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Fig. 2 Schematic Diagram of Foil Bearing Configuration 



In Refs, [l] , [5] it was shown that for © the order of unity , 

the force of each foil sector on the rotor is virtually normal since the 
tangential force component is negligible. Furthermore, the frictional 
torque and its effect on tension variations may also be neglected. With 
these approximations, the resultant forces on the rotor are (Fig. 3) 
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T ('sin o(. t + Sin c 4 eW 

) 

(2a) 





-Wi 

11 

u!* 

T k ( s ' n + Sl " ‘’U 

) 

(2b) 

The subscripts i 

and e denote in the above 

" inlet" and "exit" res- 



pectively. 

The foil length is considered to be comprised of (a) two nearly 
straight sections, (b) of a circular arc wrapping the rotor, and (c) of 
two circular arcs wrapping the guides. In addition, correction terms 
<± [k and <J e k are added to account for slack due to foil flexure (see 
Appendix A) . Mathematically 

4 * 1 k * + *Sr. )K u-AO + +<sfk + <U 

r r < 3) 

In Eq. (3) dr 0 , drg^ represent the thermal* increments in r D and r g ^ 

relative to the reference state. 

For evaluation of the thermal expansion of the structure and the 
elongation of the foil, the following thermal model will be assumed. 

The journal will be lumped at one temperature T r . Each shield-damper 
will be considered to be at a single temperature . The free length 

of the foil will be assumed to be at a temperature > whereas the 


* Ar 0 may also include the centrifugal growth of the rotor. 


6 






the portion of the foil which is in contact with the guides wiri ue assumed 
to be at the temperature 'tgy. . The excess of these temperatures over nr 0 , 
the temperature at the reference state, will be denoted respectively by 

• With t ^ ie a b° ve assumptions, the elongation of the 

foil is : 



The following geometric relations complete the formulation: 

(r.+l,** <■&!') si* *(. + <5a> 

kE^£k&k+X.L-(r o +fi*+,flZ+% Sin. “W 

+ * y„ *( r.*h k + Sf 0 * Cts et £ it - Lit sm <*;t (5c) 


i / . (5d) 

4* + - i^o * + ■ S’n + -f fj cos * e If - Kk s,n *e4 

To summarize, given the instantaneous journal position x,y and the 
temperature distribution (dictating cTa^, cTb c j c , ) and given the 

clearance h£. , the above set of 2 +(3x8) equations can be solved for the 

unknowns F^,Fy and Xj^ , yj^, T , o^» ik * ^"ek 7 ^ ik 7 ^ek 7 ^ k 
where k = 1 , 2,3. 
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3.0 SIMPLIFICATIONS 


Considerable simplification without loss of generality is obtained 
if consideration is restricted to journal displacements and thermal expan- 
sions such that: 


x fc » 3 >k > Sr ° , , Kt, , 'U ■« n. 

toil. > •* © 


Neglecting products of perturbations from the reference state, the per- 
turbation in foil length can be found from Eq. (3), 


-e(ln|-<r5\ro)© +2S^u^ 


The unknown differentials in Eq. (6) can be found by taking differentials 
of Eqs.(5) . 




cos — 

2 . 
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Using EqS 4 (7), one can deduce the relations 

"Kk * ^ek = 

4-^.6-^ (8a) 


^ct COS | -^(3k + &a cO si> '| - (n, + ^)(Hu + H 


Substituting Eq. (8) into Eq. (6), and Eq. (6) into Eq. (4) yields: 

T k -T 0 - 1 ^ Si* f + * S<T lt ♦ - M xU 

11 ■ 53 — ■ ■ ■■ ■■■ ■— i — — ■ — — — — ^ 

EJ 4ff 


Eld 

where 




■4ff- 


k 4 , - r 0 H* 


H k = 0.643 + 0.286 


i^u 1 


V r - 


+ 1-^05 


Tk/r. 


= ft k co S | - | 


*5 h ^ - Sr t 




no) 


(11) 

Ei3 

T l 

(12a) 

- 0.185 

1 k 

p.q 

(12b) 


C8] 

2 Jo. / — 

ck x 

(13a) 



(13b) 
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The above equations deserve a few comments. The variable J . r 

eff 

given by Eq. (11) describes the effective foil length, taking into account 
the fact that extension of part of the foil is restrained due to frictional 
contact with the guides. Equation (12b) gives an approximate curve fit 
to account for the effects of inertia and compressibility of the fluid. [8] 
Finally, Eqs. (13) state the assumption that the shield-dampers and their 
supporting structure are at sensibly uniform temperatures and the fact that 
their coefficients of thermal expansion are equal. 

Using simplifications consistent with the derivation of the fore- 
going equations, the bearing forces on the rotor can be reduced from Eq.(2) 
to the form: 


F, - \[3 six. | 

F y - 2 six, | (T,-^4-) 


(14a) 

(14b) 


Since axial symmetry is assumed for the present bearing, the sub- 
scripts k may be omitted in Eqs. (9) - (13). An iterative solution for T 
and h* is incorporated in the program given in Appendix B. 
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4 . 0 STIFFNESS DETERMINATION 


The evaluation of the equivalent spring constants of the foil- 
bearing system, requires differentiation of Eqs. (14) , for which the deriv- 
atives are needed. The evaluation of these expressions 

9 * ) 9 ^ 

from Eq.(9) necessitates, in turn, finding the magnitudes of 


For this purpose, use will be made of the following functional forms de- 
rived from Eqs. (5) and Eq. (A26) of Appendix A respectively: 


h- U*k . UL > r -) 


(15a) 



(15b) 


with similar expressions for bg^ . For simplicity, the dependence 
upon additional parameters which are not affected by rotor displacement 
has not been shown in Eqs. (15). It should also be noted that the shown 
dependence upon r stems from the fact that for the purpose considered in 
this section the ra8ius extends to the foil and includes h*, which is an 
implicit function of x and y . 
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Differentiation, therefore, gives: 



Similar expressions are obtained byreplacing i by e and x^ by y^ . 
Differentiating Eq. (9), with respect to x^, y^ and collecting terms 
gives respectively; 



For the particular case of no static load, the equilibrium 
position is x = y =0 and consequently, at this position 


T.-TV-T* = T 

(19a) 


(19 b) 

b ;tc = Hsf b 

(19c) 


(19d) 



(19e) 


— 1 — 9 

(19 f) 


Xk o>Cl 

^b Ltc 

(I9g) 


It follows from Eqs. (17)and(18) that under these conditions 
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This allows evaluation of 


PT W 

form: 

3 < 

> sy 


3T. 

Sx k 

> 

ii 

<V 

3* * 


Thus, 
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— f\ c,os» ^ 


by means of formulae of the 






' Sl * ty “ 

O 

’ SI* ^ - 

§ 
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x 

COS. f y = 

\ 

■ C O ^ s 

l 

"£ 


_ \ 
~ X 




(21) 


(22a) 


(22b) 


It follows from Eq. (14) that the stiffness matrix of the bearing at the point 
O is composed of the terms: 




k 

k. 


rt 


k 


v - 


25 = » s,„ ® 

3* J - 

R 

(23a) 

^-3 sw-S- 
*- 

(\ 

(23 b) 

.2£. o 





(23c) 

2fa = o 

3x 


(23d) 
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5 . 0 THERMAL ANALYSIS 


5.1 Introductory Remarks 


It was shown how the operating tension, the fluid film thickness 
and the bearing stiffness could be determined, provided that lumped temper- 
atures of the foil, the shield-damper-guide and the rotor were known. The 
present section deals with the determination of foil temperature and gives 
formulas for the evaluation of temperature-dependent fluid properties. 

In principle , one could construct a mathematical model of the 
entire machine by judiciously lumping parts of it and simulating mathemati- 
cally the variation of temperatures of the lump-regions as functions of time. 
This approach, however, would necessitate a very extensive effort. A more 
tractable analysis was to prescribe upper and lower bounds on the tempera- 
tures T g, T f on the basis of available information and to compute the cor- 
responding bounds on h* k, T, T ^. The results will be limited to the axi- 
symmetrical case. 


5 , 2 Determination of Foil Temperature 


The steady state, olanar model shown in Fig. 4 is considered. 

It is presumed that Couette flow prevails between the foil and the rotor. 
The foil-temperature variations are neglected and the temperature distri- 
bution in the rather stagnant gas between the foil and the "shield-dampers" 
is taken to be linear. With the assumption of Couette flow, the tempera- 
ture distribution across the fluid film is [9] : 


-% j- a u ‘ a (. z ) 
■ k* + k*[ fi'/ 



where \ a is the thermal conductivity of the fluid. Assuming a quasi-steady 
state (or, equivalently, neglecting the thermal inertia of the foil) the net 
heat flow into the foil is zero so that 



zl? 


Tt-'h 
\\^ - - d 


= o 


(25) 
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Fig. 4 Schematic Diagram of Temperature and Velocity 
Fields Between Rotor and Shield- Damper 



It follows that 


X 

f 


/f.U 1 h s ~lv~ d ^Iri 

Tk -\d — + * + l^r 


(26) 


where the spacing h g is a function of the thermal expansion of the rotor 
and the shield-damper , expressed by the equation: 


( % -\) u a i q - 11 - ("fr- -'c .) c 


The spacing between the rotor and the shield-damper in the reference 
state is denoted in Eq. (27) by h go . 

5.3 Variation of Fluid Properties 

The properties of the lubricating fluid are evaluated in this re- 
port at the temperature 't’a = ( Xf + T r )/2 r and at the pressure p a 
of the fluid entering the bearing. The fluid density is found from the 
perfect gas equation: 


17 


( 28 ) 


* 'B'f.fak.) 


The viscosity |i a and thennal conductivity k a are found from the 
Sutherland formulae [10, ll] 




+ C /< 



(29) 




f ^b fftU) + CA 

T 

\ ^ + C A / 

\ ^b(*kO 



(30) 


In Eqs. (29) and (30) the base viscosity |i. b and conductivity X. b 
are given at the temperature (using the absolute scale) and the 
constants C , are known from results of curve fitting over a range 
of temperatures. 
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TEMPERATURE t (°K) STIFFNESS PER UNIT BEARING WIDTH TENSION PER UNIT BEARING WIDTH FILM THICKNESS h*(IO' 8 x m) 

k OO “x N/ln/m) T(IO s xN/m) 













Geometrical Parameters 


r Q = 0.8750 in (0.02222 m) b = 0 . 1 493 in (0 . 003 79 m) 

r g = 0.1225 in (0.00311 m) d g = 0 . 1 250 in (0 . 003 1 7 m) 

r c = 1.0096 in (0.02564 m) 0 = 88.8° 

h = 0.0053 in (0.00013 m) 0 = 133.5° 


Elastic Characteristics 

T q = 2’.0 lbf/in (350. N/m) f = 0.15 

E = 30x 10 6 psi (20.7xl0 10 N/m 2 ) \> = 0.3 


Fluid and Thermal Characteristics 


Working Fluid - Air 

p a = 14 . 7 ps i (10 . 1 x 10 4 N/m 2 ) 

c* g = 5 . 6 x 10“ 6 °F“ 1 (1 .01 x 10” ^ °K _1 ) 


M = 28.97 



= 70°F (294. 3°K) 

7 . 3 x 10“ 6 °F -1 ( 1 . 31 xlO- 5 o K -l) 


Viscos ity_Parameters_(Eq. 29) 

(j. b = 0. 262 x 10~ 8 lbf • sec/in 2 (1 . 806 x 10" 5 N • s/m 2 ) 
% b = 68°F (293°K ) 

= 216 °F (1 20°K) 


ThermaJ_ (^onductivUy_ Parameters_(Eq . 30) 

X b = 0.0149 BTU/(hr- ft °F)(1. 706 J/(m.s°K)) 
= 68°F (293°K ) 

C x = 225°F (1 25°K) 


Fig. 5 Estimated Range of h* , T, k, Tf Pertinent to operation with 
Air (Acceleration, Thermal Transient at Nominal Speed, 
and Coast Down for One Combination of Foil Thickness d 
and Foil Thermal Expansion Coefficient oCf. ) 
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FOIL THICKNESS d- 0.7 mil (17.8 * IO*®m) 

COEF. OF THERMAL EXPANSION a, ■ 9 6 * IO‘®l/*F 



100 200 300 400 500 

(d> 


Geometrical Parameters 


r o = 

0.8750 in (0.02222 m) 

b 

= 0.1493 

in (0.00379 

m) 

r g = 

0.1225 in (0.00311 m) 

d g 

= 0.1250 

in (0.00317 

m) 

<c = 

1.0096 in (0 .02564 m) 

0 

= 88.8° 



II 

O 
i Cn 
-C 

0.0053 in (0.00013 in) 

®g 

= 133.5° 




Elastic Characteristics 

T 0 = 2.0 lbf/in (350. N/m) f = 0.15 

E = 30x 10 6 psi (20.7x 10 10 N/m 2 ) V = 0.3 


Fluid and Thermal Characteristics 

Working Fluid - Air M = 28.97 

p a = 14.7 psi (10.1 xlO 4 N/m 2 ) ^room = 70°F (294 . 3°K) 

<* g = 5 . 6 x 10" 6 °F" 1 (1 .01 x 10“ ^ °K -1 ) et f = 7 . 3 x 10" 6 °F _1 (1. 31 x 10” 5 °K~ 


Viscqsity_Parameters_(Eq. 29) 

(±5 = 0. 262 x 10 -8 lbf- sec/in 2 (1 . 806 x 10“^ N • s/m 2 ) 
^ b = 68°F (293 °K ) 

= 216 °F (1 20°K) 


Thermal^ Conduct_ivity_Parameters_(Eq. 30) 

\ b = 0.0149 BTU/(hr-ft °F)(1. 706 J/(m.s°K)) 
% b = 68°F (293°K) 

C x = 225°F (1 25°K ) 


Fig. 6 Estimated Range of h*, T, k, Tf Pertinent to Operation 
with Air (Thermal Transient at Nominal Speed, for Four 
Combinations of Foil Thickness d and Foil Thermal 
Expansion Coefficient o(f • ) 
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TEMPERATURE T CK) STIFFNESS PER UNIT 8EARING WIDTH TENSION PER UNIT BEARING WIDTH FILM THICKNESS h* (I0' 8 x m) 

k (IO B x N/Wml T (I0 3 x N/m) 

















b = 0.1493 in (0.00379 m) 

d g = 0.1250 in (0.00317 m) 

© -- 88 . 8 ° 

© n = 133.5° 


Elastic Characteristics 

T q = 2.0 lbf/in (350. N/m) f = 0.15 

E = 30x 10 6 psi (20.7 xlO 10 N/m 2 ) V = 0.3 


Fluid and Thermal Characteristics 

Working Fluid - Xe-He Mixture M = 83.8 

p a = 25. 1 psi (17.3xl0 4 N/m 2 ) t r oom = 70°F (294 . 3°K) 

<X g = 5.6xl0 -6 OF -1 (l.OlxlO -5 °K -1 ) <*. = 7.3 xlO-6 op-1 ( 1>31xl0 -5 o K -l) 


Videos Uy_P^rajrieters_(Eq. 29) 

PJ-, = 0. 2475 x 10~® lbf - sec/in 2 (1.706 xl0 - ^ N*s/m 2 ) 
T b = - 1 00°F ( 200°K ) 

C M = 34 6 . 74 °F (192. 66°K) 

r 


Thermaj_ (^on_d^cUvj_ty_Param_eters_(Eq . 30) 

\ b = 0.0110 BTU/(hr-ft°F) (1 .056 j/( m . s °K) ) 
t b = - 1 00°F ( 200°K ) 

C, = 168.1 6°F (93 . 42°K) 

A. 


Fig. 7 Estimated Range of h* , T, k, Tj Pertinent to Operation 
with Xe-He Mixture (Acceleration, Thermal Transient at 
Nominal Speed, and Coast Down for One Combination of 
Foil Thickness d and Foil Thermal Expansion Coefficient o£f.) 


Geometrical Parameters 

r G = 0.8750 in (0.02222 m) 
r g = 0.1225 in (0.00311 m) 

r c = 1.0096 in (0.02564 m) 

h o = 0.0053 in (0.00013 m) 
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ESTIMATE -(C) 




















300 350 400 450 500 



Geometrical Parameters 


r o “ 

0.8750 in (0.02222 m) 

b 

0.1493 in (0.00379' m) 

r g = 

0.1225 in (0.00311 m) 

dg = 

0.1250 in (0.00317 m) 

r c = 

1.0096 in (0.02564 m) 

0 

88.8° 

zr 

O 

11 

0.0053 in (0.00013 m) 

©g = 

133.5° 

Elastic 
T = 

Characteristics 
2. 0 Ibf/in (350. N/m) 

f = 

0.15 

E 

30x 10 6 psi (20.7xl0 10 

N/m 2 ) V - 

0.3 


Fluid and Thermal Characteristics 


Working Fluid - Xe-He Mixture 
p a = 25.1 psi (17.3xl0 4 N/m 2 ) 

<X g = 5 . 6 x 10“ 6 °F _1 (1.01 x 10“ 5 °K _1 ) 


M = 83.8 

'troom = 70°F (294. 3°K) 

tX. = 7.3x10-6 o F -l (i. 3ixl0- 5 °K _1 ) 


ViscgsUyJ 3 arameters_(Eq. 29) 

Hb = 0 . 2475 x 10 -8 lbf - sec/in 2 (1.706 x 10" 5 N ■ s/m 2 ) 
r b = - 100°F ( 200°K ) 

C = 346 . 74 °F (192 . 66°K) 



Thermal_ <^ond(^cUvUy_Param^ters_ (Eq . 30) 

X b = 0.0110 BTU/(hr-ft°F) (1.056 j/(m-s°K) ) 
T b = - 1 00°F ( 200°K ) 

C x = 168.1 6°F (93 . 42°K ) 


Fig. 8 Estimated Range of h*, T, k, 'tf Pertinent to Operation 
with Xe-He Mixture (Thermal Transient at Nominal Speed 
for Four Combinations of Foil Thickness d and Foil 
Thermal Expansion Coefficient e £{- ) 
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6.0 RESULTS 


A computer program has been developed to give h*, k, T, 'tf 
(referred to below as the "output parameters") for a given set of input 
parameters. Scanning capabilities for certain input parameters are built 
into the program. These capabilities are explained by comments in the 
program listing (Appendix B) . 

6 . 1 Method of Presentation 

Values of the output parameters are shown for a range of con- 
ditions in Figs. 5-8. Estimates pertinent to operation with atmospheric air 
as the working tluid are presented in Figs. 5 and 6. In analogous fashion, 
estimates pertinent to operation with a mixture of xenon and helium (83.8 
molecular weight) are given in Figs. 7 and 8. For the particular combina- 
tion of foil thickness and thermal expansion coefficients used in the design. 
Figs. 5 and 7 show estimates pertinent to (a) start-up with r<apid accelera- 
tion, (b) slow thermal transients which follow, with the machine at rated 
speed, and (c) coast-down. Figures 6 and 8 show a comparison of estimates 
pertinent to thermal transients for four combinations of foil thickness and 
coefficient of thermal expansion. (The first column in each figure pertains 
to the present design.) 

(a) Results pertinent to rapid acceleration (a few seconds) 

In order to study the bearing characteristics during start-up 

the output parameters are found as a function of operating speed 

with T and T assumed to remain substantially constant. These 
r g 

assumptions are justified by the rapidity of the process. The tem- 
perature of the low-thermal-inertia foil may range between its 
initial temperature (for infinitely fast acceleration) and the temper- 
ature predicted by Eq. (26) for sufficiently slow acceleration) . The 
results are shown in Fig. 5a (air) and in Fig. 7a (Xe-He) . At 600 
rps , the stiffness predicted on the basis of the above assumptions 
is in the neighborhood of 20,000 lb/in/in (14 x 10^N/m^) and the 
film thickness is approximately 250 microinches for air (6 microns) 
and 350 microinch for Xe-He (9 microns) . 
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(b) Results pertinent to slow thermal transients at rated speed 

Following start-up, while running at rated speed, the system 
warms up relatively slowly to its equilibrium temperatures. 

The shield-damper-guide temperature follows that of the rotor. 

A reasonable range of rotor and guide temperatures (based on 
available thermal maps [12] of the pivoted-shoe bearing 
version of the machine) , is used to predict the range of the 
output parameters. It is assumed in this case that the foil 
temperature changes quasi-statically according to Eq. (26). 

All results are plotted versus the rotor temperature T r which 
may be regarded as a distorted time scale. By comparing Figs. 

6a and 6b, 6c and 6d, 8a and 8b, 8c and 8d, one may observe 
that a reduction of foil thickness from of 1.0 to 0.7 mils (25 
to 18 microns) causes little change in bearing stiffness and gap. 
On the other hand, comparison of Figs. 6a and 6c, 6b and 6d, 

8a and 8c. 8b and 8d indicates that changes in foil material 
L thermal expansion coefficient 6 . 3 x 10”® o -1 (11.3x10”® 
o^-l) to 9.6 x 10”® op-1 (17.3 x 10 “6 o^-l) J may cause signifi- 
cant change in bearing behavior.* If the coefficient of thermal ex- 
pansion of the foil is too low relative to that of the journal and of 
the foil support, the gap may be reduced excessively. Conversely, 
too high a value of the coefficient may produce a large clearance 
and an overly loose bearing. In the assumed ranges of temper- 
atures, Figs. 6a and 8a confirm the correct selection of materials 
and other design parameters and predict favorable operation in 
terms of adequate st'ffness and safe film thickness.** 


* The curves in Figs. 8c and 8d were not carried to the upper temperature 
range due to convergence difficulties. It was not deemed worthwhile to 
make special program provisions for this lunge, sinue in this extremely 
relaxed foil state the accuracy of Eq. (29b) is marginal. 

** Coefficients of thermal expansion of foil and foil supports were matched on 
the basis of a simplified calculation carried out by L. Licht prior to the 
present analysis. Metals of compatible thermal expansion and suitable 
elastic properties were then selected accordingly. The manufacture of foil 
bearings for the BRU has been completed at the time of drafting of this 
report. 
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(c) Results pertinent to coast-down 

Coast-down is moderately rapid. Therefore, it is plausible to 
assume that during this process, the journal and the foil sup- 
port do not experience appreciable changes in temperature. On 
the other hand, the thin foil is expected to follow quasi-statical- 
ly the temperatures dictated by its more massive neighbors and 
by the heat generated in the film. Thus, the foil temperature is 
estimated on the basis of Eq. (26) for a range of rotor and foil 
support temperatures. The results are shown in Fig. 5c (Air) 
and Fig. 7c (Xe-He) . The results indicate an increase in stiff- 
ness duriny coast-qown, corresponding to a 10% increase in 
resonant frequency in comparison with start-up. 
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7.0 CONCLUSIONS 


The parameter values selected in the design of the BRU foil 
journal-bearings have been tested by the computer program developed in 
this report. For the temperature ranges assumed for the journal, foil and 
foil supports, the design has been confirmed with regard to bearing stiff- 
ness and film thickness. Since foils can be easily replaced, the foil 
material is a sensitive variable for controlling the bearing characteristics, 
which remains at our disposal when other parameters cannot be altered. 
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APPENDIX A 


EFFECT OF FOIL FLEXURE 

In this appendix, a model is presented for the calculation of 
slack resulting from the nonvanishing foil stiffness. 

Cons Lder two parallel cylinders of radii r Q and r g res- 
pectively (Fig. Al). The center to center distance of the two cylinders 
is r c . The corresponding distance b is given by the equation 



Fig. Al Schematic View of a Foil Wrapped 
Around Two Cylinders 
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Let the foil be wrapped around the cylinders at a tension T The slack 
Cj is defined as the excess of the length AqA^A 2 A 3 of a finite-stiffness 
foil over the length AqA^A;^ of a perfectly-flexible foil. In order to 
calculate (j, the contour of the foil must be found. 

Let a coordinate s be measured along A A,A„A„ with A 

o 1 2 3 o 

as the origin. At any point s, the distance w to the corresponding po- 
sition of the perfectly-flexible foil is measured normal to the foil of finite 
stiffness. The local radius of curvature of the foil is R. R is taken posi- 
tive when the curve is concave down (Fig. Al) . Since there are no normal 
loads on the foil the differential equation for the curvature is 



If the slope of the perfectly flexible foil does not differ much 
from that of the finite-stiffness foil, so that 



it follows that 

I ~ _1 ct*ur 


(A3) 


(A4) 


where R^ is the local radius of curvature of the perfectly flexible foil 
at s . 

It is convenient to use the following dimensionless notation. 


a 


D 


(A5) 


W 


u r 

r 0 e x 


(A6) 
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s 


(A7) 


The problem is then formulated by the following differential equations: 


where 



r 


“31 V 




V- p- 

3 




(A 8 ) 


(A9) 


and ^j ,%2 ,^3 correspond to A^Az,^ in Fig. Al, and 



* 


The boundary conditions are: 
at -O 


(A10) 



(All) 


W =0 

Jul 



(A12) 


(A13) 
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at 





15 = ‘ *5 

(A14) 

£ 

li 

o 

(A15) 

dW 


»v =0 

(A16) 




An additional constraint is 



(A17) 


The above formulation contains four unknown integration constants and 
three unknown parameters * The seven conditions (All) to (A17) 

are used to determine these unknown constants. When the Eqs. (A8) through 
(A10) are integrated in the three regions^ Eq. (A9)] , and the values of W 
and dW/d£ are matched at points Aj,A 2 the solutions are found to oe: 


-A sink ^ 

U = <-A sink ^ 
-A si*iln§ 


— Cosh ^ 1 4 - A£ + I 

-Cosing ■+• + t - 

— ^ +c^ *E 


- CoS 


5 2 - 


where 



(A18) 


(A 19) 
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C - $ + /l Cosh § + SinUfc 


(A20) 


J f A sirth ^ + Cosl\ ^ — c| 


(A21) 


The values of £3 must be found from the simultaneous solution 

of the equations: 


■7^ + A(cosV.^-i) + s'.mU ^ 


(A 2 2) 


2 + ^~2Coslv^4-2A£ a C*>sh^+2^Si*U^ 

^ (A23) 

y ~ x- (A24) 


The sought excess length Cf may then be found from: 



(A25) 


(A26) 
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APPENDIX B 


PROGRAM TO EVALUATE FOIL SEARING OPERATING TENSION T, GAP H STAR AND 
C STIFFNESS CAY FUR ARBITRARY PARAMETERS. INITIAL INPUT DATA IS 
C COMPILED .IN HEJNG SUBROUTINE. DESIRED MODIFICATIONS IN PARAMETERS 
C SHOULD BE MADE THROUGH USE OF NAMEL 1 ST/ DAT A/ . CALCULATION OF T, HSTAR 
C AND CAY IS MADE IN SUBROUTINE EQN1 . AUX I L I AR Y COMPUTATION OF SLACK 
C SIGMA AND ITS DERIVATIVES IS MADE IN SUBROUTINE SLACK. 

C COMPUTATION OF FOIL TEMPERATURE AND TEMP. - DEPENDENT FLU ID— 

C VISCOSITY, CONDUCTIVITY AND DENSITY IS IN SUBROUTINE TEMPF. 

C PROGRAM OUTPUT IS IN (1) BRITISH UNITS AND (2) I NT ERNAT IQNAL < S I ) UNITS. 
C BRITISH UNITS USED IN PROGRAM ARE I NCH-L BF- SEC-DEGF OR AS 
C SPECIFIED. BASIC SI UNITS ARE M-KG-S-K. ALSO J,N. 

C E.G. VISCOSITY IN NS/ M2 ( NEWT ON*S ECOND/MET ER**2 ) 

C VARIABLES USED FOR SI HAVE I NT ENDING AND ARE PARENTHESISED (BELOW) 


C INPUT OF PARAMETERS IS IN BRITISH UNIS. DEFINITIONS: 

C RPS = REVOLUTIONS PER SECUND OF ROTOR 

C RO •= ROTOR RADIUS (PROGRAM ADDS TO THIS THIKNS/2.) ( ROT NT) 

C RG = GUIDE RADIUS (PROGRAM ADDS TO THIS THIKNS/2.) ( RGINT ) 

C RC = DISTANCE OF CENTER OF GUIDE FROM CENTER OF ROTOR (RCINT) 

C SPACE = COLD SPACING BETWEEN ROTOR AND KIDNEY SHAPED GUIDE (SPCINT) 
C TZ = INITIAL TENSION PER UNIT WIDTH (PRELOAD) (TZINT) 

C THI KNS=FO IL THICKNESS INCH (TKINT) 

C E = MODULUS OF ELASTICITY PSI (EINT) 

C ENU = POISSON RATIO OF FOIL 


C 0= BENDING RIGIDITY OF FOIL. CALCULATED INTERNALLY FROM E AND 
C THIKNS. IF 0 IS PRESCRIBED ZERO, FOIL IS ASSUMED PERFECTLY FLEXIBLE 


C FCUEF = FRICTION COEF BETWEEN FOIL AND GUIDE 
C EMOL = MOLECULAR WEIGHT OF GAS 

C PA = AMBIENT PRESSURE PSI (PAINT) 

C ALFAF = COEF. OF EXPANSION OF FOIL 1/DEGF (AFINT) 

C ALFAG = COEF. OF EXPANSION OF GUIDE 1/DEGF (AGINT) 

C ALFAR = COEF. OF EXPANSION OF ROTOR I AR I NT ) 

C TROOM = ROOM TEMPERATURE DEG F (TRMINT) 

C TGUIDE= FOIL GUIDE TEMPERATURE ( TGINT) 

C TROT GR= ROTOR TEMPERATURE (TRINT) 

C TF OIL- FOIL TEMPERATURE. I F ( I FO I L . EQ . 1 ) PROGRAM COMPUTES TFOIL(TFINT) 
C IFGIL IS ASSUMED TO BE 1 UNLESS OTHERWISE SPECIFIED 

C COND = INPUT THERMAL CONDUCTIVITY OF GAS BTU/I HR*FT*DEGF ) (CNOINT) 

C V I S= INPUT VISCOSITY OF GAS L0F*S EC/ I N**2 (SEE CONSPR) (VISINT) 


C CONSPR = LOGICAL VARIABLE. ASSUMED 'TRUE UNLESS OTHERWISE SPECIFIED 
C IF TRUE, INPUT OF VIS AND COND IS ASSUMED CONSTANT 
C IF FALSE THESE INPUTS ARE ASSUMED TO BE AT REFERENCF 
C TEMPERATURES, TRVISC AND TRCOND , RESPECTIVELY. GAS 
C TEMPERATURE ASSUMED AT AVERAGE OF TROTOR AND TFOIL. INTERNAL 
C VISCOSITY (EMU) AND AND CUNDUC T I VI TY ( CONOCT ) ARE ,THEN, CALCULATED 
C BY THE SUTHERLAND FORMULA. CVISC AND CCOND ARE REQUIRED INPUTS OF 

C SUTHERLAND CONSTANTS 

C GOMTRY= LOGICAL VARIABLE. IF TRUE, OUTPUT FORMAT FOR GEOMETRY 
C VARIATIONS. 

C T£RMU= LOGICAL VARIABLE. IF TRUE, OUTPUT FORMAT FOR TEMPERATURE 
C VARIATIONS. 

C FLUID = LOGICAL VARIABLE. IF TRUE, OUTPUT FORMAT FOR CHANGES 
C IN FLUID PROPERTIES 
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C SOLID = LOGICAL VARIABLE. IF TRUE, OUTPUT FORMAT FOR SOLIO 
C PROPERTY CHANGES 

C IF ONE SPECIFIES GOMTRY =.TRUE. AND AT LEAST ONE CARD 
C LATER, ONE SPECIFIES GOMTRP= .FALSE. , THE COMPLETE 
C NAMELIST IS PRINTED. SAME WITH FLUIDP =. FALSE. 

C FOLLOWING FLUIDS. TRUE. SAME WITH SOLIDP AND 
C TERMOP FOLLOWING SOLID AND TERMO RESPECTIVELY. 

C SKI P=. TRUE. : LEAVES ONE LINE SPACE BEFORE OUTPUT AND SWITCHES 
C SKIP BACK TO .FALSE. 

C REAORC= .FALSE. : THETA1 (DEG) MUST BE PRESCRIBED, RC IS CALCULATED. 
C IF I SC A N=4 OR 5 READRC IS SWITCHED TO TRUE AFTER FIRST CASE 
C SCANNING OPTIONS: INITIALLY, PROGRAM ASSUMES 

C ISCAN=1 1.6. EACH NAMELIST RECORD GIVES I OUTPUT RECORD 
C I SCAN= 2 SCANS THIKNESS IN SOLID VARIATION FORMAT IN 
C (UNLESS OTHERWISE SPECIFIED) STEPS OF DTHIK=0.0002 
C TO THIKMX=0.002 

C I SC AN= 3 SCANS PRELOAD TZ IN SOLID VARIATION FORMAT IN 
C (UNLESS OTHERWISE SPEC I FI ED) STEPS OF DTZ=0. 5 TO TZMX=4. 

C I SC AN=4 SCANS RO,RC WITH CONST RADIAL SPACING IN GEOMETRY 
C VARIATION FORMAT IN STEPS (UNLESS OTHERWISE SPECIFIED) QF 
C ORO =0.25 TO R0MX=2. 

C I SC AN=5 SCANS RG,RC WITH CONST RADIAL SPACING IN GEOMETRY 
C VARIATION FORMAT IN STEPS (UNLESS OTHERWISE SPECIFIED) OF 
C 0RG=O. 05 TO RGMX=0 . 3 

C ISCAN =6 SCANS TROTOR AT RPS=CGNST, AND UNLESS OTHERWISE 
C SPECIFIED AT STEPS DTRGT=40. TO TRUTMX=270., TGIDMX=200. 

C TGUIDfc = TROTOR-C TROTMX - TGIDMX) * (TROTOR - TROQM)/ (TROTMX- TROOM) 

C TERMO OUTPUT FORMAT 

C l SC AN= 7 SCANS RPS AT CONST TEMP. DRPS=25. RPSMX=200. 

C UNLESS OTHERWISE SPECIFIED 

COMMON THIKNS, RADIUS, DG,RGUIDE,RC,DD, SPACE, BB,RG,RO , 

1 T, TZ, TT, ET, TZP,TP, EPS , EPS23, V , RPS, Y . CAY, 

2 HSTR, HSTAR , ETU, PIR2, THETA, THETA1 , TH2 , 

3 CTH2, STH2 , THETAG, IHETG1, P I , DEG , FOA , FG A , 

4 PA IE, COMP, DENSTY, PA, EMU, CONDCT CF , EMQL , C,AIE, 

A CUNSPR,TRVISC,TRCOND,CVISC, CCOND , CONI', VIS, 

5 ELTOR, RGBS , ALF , S I GBRO , S I GBR, DBDYK, DSGDB, DSBRDT , DSGDR , 
B DBDR , 

6 THERML, ALFAF, ALFAG, ALFAR, TROOM, DTF , DTG , DTR, 

3 TROTOR ,DTROT, TROTMX, TGUIDE, TGIDMX ,TFOIL, 

9 DTHIK, THlKMX,DTZ»TZMX»DRO,ROMX»DRG , RGMX ,DRPS,RPSMX» 

7 PEXT , E, FCOEF , ENU, D , SMAX 

LOGICAL TERMOP, FLU I OP , SOL I DP , GOMTRP ,CONSPR 

LOGICAL TERMO, FLUID, SOLID, GOMTRY .SKIP .READRC 

COMMON TERMOP, FLUIDP, SOL IOP, GOMTRP , ISCAN, IFOIL 

COMMON TERMO, FLUID, SOLIO, GOMTRY .SKIP, READRC 

COMMON AT ( 5 ) , AX3(5), ASIGBR(5), ABSA(5) 

COMMON AFP, AGP, ARP, AFINT, AGINT, ARINT, AFINTP, 

1 AGINTP , ARINTP, VISCP, TINT, HINTP, HINT, 

2 CYINTP, CYINT, SMINTP, SMINT, ROINT, RG I NT, RCINT, 

3PAINT, PAINTP, TRMINT, ONSINT, VISINT, VSINTP, CNDINT, 

4TKINTP, TKINT, EINT, TZINT, TFI NT , TGI NT , TR I NT , EP , £ I N TP, 

5 TCNV, CAYCNV, DNSCNV, VISCNV .CONCNV 

CALL HEDNG (1) 

1 CALL EQN 1 

CALL HtDNGI 2 ) 

GO TO 1 
END 
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SUBROUTINE HE PNG (J) 

NAMEL IST/DATA/THIKNS , OTH'IK, THIKMX,TT, ET, 

1RAOIUS* RO,DRO,ROMX,RGUIDE,RG,DRG,RGMX, RC, Y, SPACE, 

2DG, 3B,ALF,RGBB, ELTOR, 

3 THETA, THETA I ,TH2 ,CT H2 , S T H2 , THET AG, THETG 1 , 

AT , TZ , DTZ ,TZ MX , TZP , TP , 

5 P EXT , E, FCOEF , ENU » D, ALFAF, ALFAG, ALFAR * 

6PAIE ,COMP»OENSTY , PA, VIS » EMOL , COND.EMU, CONDCT, TR V I SC, TRCOND , 
7RPS , DRPS, RPSMX, V, ETU, PIR2, CVISC.CCONQ, 

8 TF 01 L , TGU I OE , TROTOR , OT ROT , T ROTMX , TG I DMX , TROOM , THERML , 

9TERMC1 ,FLUl 0 , SOL If.) , GOMTRY , TER MOP ,FLUIUP,SOLIQP, GOMTRP , CONS PR, 
AISCAN, IFOIL ,SKIP, READRC 

NAMEL I ST/SI /TKINT »ROINT ,RGINT ,RCINT ,T I NT , TZINT, SPCINT, 

1 E I NT , AFINT, AG I NT, AR I NT, DNSINT, PAINT, CNQINT , 

2VISIMT, TF I NT , TG I NT , TR I NT , TRMINT 
COMMON THIKNS, RADIUS, DG, RGU I DE ,RC , 00, SPACE, BB,RG,RO , 

1 T, TZ, TT, ET, T Z P , T P » EPS , EP S23 , V , RP S ,Y , CAY, 

2 HSTR, HS FAR, ETU, PIR2, THETA, THETA1, TH2, 

3 CTH2, STH2 , THETAG, THE TGI > P I , DEG, FO A, FG A, 

4 P A I E , COMP, UENSTY, PA, EMU, CONDCT, CF, EMOL, C,AIE, 

A CUNSPR , TR VI SC , TRCOND , CV ISC, CCOND, CONO,VIS, 

5 ELTOR, RGBB , ALF , S I GBRO , S I GBR , DBDYK, DSGDB, DS8RDT,DSGDR , 
B DBDR, 

6 THERML, ALF AF , ALFAG, ALFAR, TROOM, DTF, DTG, DTR , 

8 TROTOR , DT ROT , TROTMX, TGUIDE , TGIDMX ,TFOIL, 

9 UTHIK, THIKMX.DTZ.TZMX , DRO , ROMX , DRG , RGMX, DRPS, RPSMX, 

7 PEXT, E, FCOEF, ENU, D ,SMAX 

LOGICAL TERMQP, FLU I OP , SOL I DP , GOMTRP ,CONSPR 

LOGICAL TERMO, FLUID, SOLID, GOMTRY ,SKIP , READRC 

COMMON TER MU PjFLUIDP, SO LlDP, GOMTRP ,ISCAN, IFOIL 

COMMON TERMO, FLUID, SOLID, GOMTRY , SKIP, READRC 

COMMON AT ( 5 ) , AX315), AS I GBR (5 ) , ABS A( 5 ) 

COMMON AFP, AGP, ARP, AFINT, AGINT, ARINT, AFINTP, 

1 AGINTP, ARINTP, VISCP, TINT, HINTP, HINT, 

2 CYINTP, CYINT, SMINTP, SMINT, ROINT, RGINT, RC I NT , 

3PAINT, PAINTP, TRMINT, DNS I NT , VISINT, VSINTP, CNDINT, 

4TKINTP, TKINT, EINT, TZINT, TFINT, TGINT, TR INT , EP, E INTP , 

5 TCN V , CAYCNV, DNSCNV, VISCNV ,CONCNV 

GO TO ( l ,2 ) , J 

I CONTINUE 

T C N V = 4.4482/2. 54E-2 

CAYCNV= (4.4482/2.54E-2**2 )*l.E-9 

DNSCN V= 386.* 2.768E4 

VISCNV- 144.* 47. 88*1. E4 

CF = 778.16/3600. 

CONCN V— 144. *5. 188 7E 2/778. 16 
C=100. 

AI E=0. 2 

PI= 3 . 1 4 1 59 26 536 
DEG= 180. /PI 
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CCCCCCCCCCCCCCCCCCC INITIALIZATION OF DATA CCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C ASSUMED INITIAL VALUES. MAY BE MODIFIED BY STATEMENT IIO 

THIKNS = O.OOl 
D TH I K=0 . 0002 
THIKMX^O. 002 
EMOL = 83.8 

RPS = 800. 

ORPS= 25. 

RP SMX=2 00 . 

PA = 25.1 
ALFA F= 6.3 E-6 
ALFAG = 5 . 6 E — 6 
ALFAR = 7.3E-6 
TROOM = 70. 

TF 01 L=470. 

TGUI0E=330. 

TGI 0MX=200. 

T R0T0R=400. 

DTR0T=40 . 

TR0TMX=270. 

E = 30. E6 
ENU = 0.3 
FCOE F =0 . 1 5 
VIS =5.lE-9 
CVISC= 346.792 
CCOND = 168.176 
T RCQND= -100. 

TR VI SC=— 100 . 

CON SPR= .TRUE. 

COND = 0.0123 
RO = 0. 875 
DR0=0 .25 
R0MX=2. 

RG= 0.1225 
DR G= 0 .0 5 
RGMX=0. 3 
RC = 1.0096 
SPACE=0. 00525 
TZ = 2.0 
0TZ = 0 .5 
TZ MX=4. 

I FOIL = 1 
ISCAN=1 

GOMTRY=. FALSE. 

GOMTRP=. FALSE. 

TERMO=. FALSE. 

TERMOP=. FALSE. 

FLU I u= .FALSE. 

FLUIOP=. FALSE. 

SOL 1 0= .FALSE. 

SOLIGP=. FALSE. 

SKIP=. FALSE. 

REA0RC=. TRUE . 

CCCCCCCCCCCCCCCCCCC END OF DATA IN I T IAL IZAT ION CCCCCCCCCCCCCCCCCCCCCCCCC 
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1= 2.*TZ 
0=0.003 
F0 A = 2. /DEG 
FGA = 2 . /DEG 
HSTR = 0.0005 
OT F= TFOIL-TROOM 

2 GO TO ( 110,120,130,140,150,160,170 ),1SCAN 

120 THIKNS=THIKNS + DTHIK 

IF(THIKNS.GT.THIKMX) GO TO 110 
GO TO 118 
130 TZ= TZ+DTZ 

IF(TZ.GT.TZMX) GO TO 110 
GO TO 118 
140 R0=R0+DR0 
RC = RG+ DRO 

IF (RO.GT.ROMX) GO TO 110 
GO TO 118 
150 RG=RG+ORG 
RC = RC +ORG 
READRC=.TRUE. 

IF (RG.GT.RGMX) GO TO 110 
GO TO 118 

160 TRQTGR= TROTOR + DTROT 

TGUIDE=TROTOR-(TRGTMX-TGIDMX)+( TROTQR-TRUOM ) /( TROTMX— TROOM) 

I F ( TROT OR. GT. TROTMX) GO TO 110 
GO TO 118 
170 KPS=RPS+ DRPS 

IF (RPS.GT.RP SMX ) GO TO 110 
GO TO 118 

110 READ ( 5 , OAT A , EN0=40 1 ) 

GO TO (118,123,128,148,148,168,168,168 ),ISCAN 

128 SOL I 0=. T RUE . 

GO TO 118 
148 GO,MTRY=.TRHE, 

GO TU 118 
168 T ERMO= .TRUE . 

118 DTG=TGUIOE-TROOM 

RADIUS = RO + THIKNS/2. 

RGUIDE = RG + THIKNS/2. 

DG= RG+ 0.0025 
RGPRO=RG+RO 
I F ( REAORC) GO TO 8 
THETA-THETA 1 /DEG 

RC=SURT< DG**2 +<RGPR0)**2 -2 . * S IN ( P I / 3. - THE TA /2 . )* ( RGPR 0) *DG) / 
1 COSIPI/3. - THETA/2.) 

GO TO 9 

3 THETA =2 . *P I / 3 . - 2. * ( ARCOSU RGPRO ) /RC ) ♦ 

1 ARSIN < OG/RC)) 

THET Al= THETA *DEG 
9 TH2 =THETA /2.0 
STH2 = SIN( TH2 ) 

CTH2 = COS ( T H2 ) 
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AFP=ALFAF*1 .Eb 

AGP=ALFAG*1.E6 

ARP =ALFAR* 1 . £6 

AF I NT P= ALF AF*1 . 8 £6 

AGINTP= ALFAG*1 . 8E6 

ARINTP=ALFAR*1.8E6 

AF I NT =AF I NTP*1 • £—6 

AGINT= AGINTP<=1.£— 6 

ARINT= ARINTP*l.E-6 

R0INT=R0*2. 54E-2 

RGINT=RG*2.54E-2 

RCINT=RC*2.54E-2 

SPCINT = SPACE *2. 54E— 2 

PAINTP=PA*6 .8947 £— 1 

C PAINT=PA*6.8947E3* (l.E-4 FOR PRINTING) 

PAINT=PAINTP* I. £4 

T RM I N T= ( T ROOM + 459.67J/1.8 

TKI NTP= TH I KNS’ it 2 . 54E4 

TKINT = TK I NTP* 1 . £—6 

E I NTP= E* 6 .8947 £—7 

EINT=EINTP*l.E10 

EP= E*1 . E— 6 

TZ I NT= TZ*TCNV 

TGINT=(TGUIDE*459.67)/1. 8 

TRINT=(T ROT OR +4 59. 67)/ 1.8 

OTR=TROTOR— TROOM 

CONDCT= CONDUCE 

TT = THI KNS*1000. 

E T= E*TH I KNS 
TINT = TCNV*T 
T P= T/ET 

IF (RC.GT.(RADIUS*RGUIDE*0.005) ) GO TO 4 
PRINT 10, RC , RO , RG 

10 FORMAT ( 26H GUIDE TOO CLOSE TO ROTOR / 10X , 2HRC , 18X, 2HRC, 18X, 
1 2HRG/ 3E20.5) 

GO TO 2 

4 BB = SORT ( RC**2- ( RADIUS* RGU I DE ) **2 ) 

DBOYK= - STH2 + CTH2* (RADIUS * RGUIDE )/BB 

DBDR=-( RADIUS *RGU IDE )/BB 

RGBS = RGUIDE/BB 

ALF = 2 . *B 8 * RADIUS *THETA 

C ALF = FREE-FOIL LENGTH 

ELT0R=ALF/RA01US 
PIR2=PI*RA0IUS*2.0 
V= RPS^PI R2 

THETA G=AR CO S(OG/( 2 . *RGU IDE)) +ARCOS( DG/RC )-ARCOS( ( RG+RO ) /RC ) 
THE TGI = THETAG*DEG 

C Y= (ACTUAL Y) /RADIUS SGBR = SIGMA/RADIUS 

Y =0. 

I F ( D. EO. 0. ) GO TO 5 

D= E*TH I KNS **3/( 12. *( l.-ENU**2) ) 

11 = 1 
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CALL SLACK < TZ,S1GBR0,RADIUS,RGU IDE, RC, 1 1 > 

I 1=1 

CALL SLACK ( T l , S I G6R0 , R AD I US , RGU I DE ,RC , 1 I ) 

5 00= SPACE - 0TR*ALFAR*R0 * QTG*ALFAG* (RC-RG) - THIKNS 

CALL TEMPF 

COMP = PA ^RADIUS/ ET 

PE X T =2. *(RGU1 DE/RADIUS )* ( 1 . — EXP (~FCOEF*THETAG) )/FCOEF 
T ZP= TZ/ET 

I F ( GGMTRY. AND.. NOT. GOMTRP) GO TO 200 
IF ( FLUID. AND.. NOT. FLUIOP) GO TO 210 
IF ( SOLID. AND.. MOT. SOL I DP) GO TO 220 
IF ( TERMO. AND.. NUT. TERMGP) GO T 0 230 
GO TO 300 

200 GOMTRP=. TRUE. 

WRITE (6,400) 

WRITE (6, DATA) 

WRITE (6, SI) 

FLUID = .FALSE. 

FLU IDP=. FALSE. 

SOLID = .FALSE. 

SOLIDP=. FALSE. 

TERMOP=. FALSE. 

TERMO = -FALSE. 

PRINT 201 

201 FORMAT ( SOX ,35HVAR I AT ION OF GEOMETRY // 

1 60H RPS T H* GAP K SMAX , 

2 60H RO RG RC THETA THETAG/ 

3 6 OH LB/IN N/M MINCH E-6*M LB/IN2 E9*N/M2 PSI E8*, 

4 60HN/M2 IN M IN M IN M DEG DEG ) 

GO TO 300 

210 GOMTRP=. FALSE. 

WRITE (6,400) 

WRITE (6, DAT A) 

WRITE (6, SI) 

GOMT RY= .FALSE. 

FLUI DP=. TRUE. 

SULIDP=. FALSE. 

SOLID = .FALSE. 

TERMOP=. FALSE. 

TERMO = .FALSE. 

PRINT 211 

211 FORMAT ( 50X, 35HVAR I ATI ON OF FLUID PROPERTIES // 

1 60 H RPS T H* GAP K PA , 


2 

60H 

TROOM 

M 

DENSITY 

VISCOSITY 

CON 

3 

10HDUC T IV1 TY 

/ 




4 

60 H 

LB/IN 

N/M 

MINCH E-6*M 

LB/IN2 E9*N/M2 PSI 

E4* 


5 60HN/M2 F K LBF*S2/IN4 KG/M3 E-8*REYN E-4*NS/M2 B/H*F, 

6 10HT*F J/MSK ) 

GO TO 300 

220 GOMIRP=. FALSE. 

WRITE (6,400) 

WRITE (6, DATA) 
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WRI TE ( 6 , SI) 

GGMTRY= . FALSE. 

FLUIDP=. FALSE. 

FLUID = -FALSE. 

SOL IDP=. TRUE. 

TERMOP=. FALSE. 

TER MO = -FALSE. 

PRINT 221 

221 FORMAT ( 50X , 3 5HV AR I AT I ON OF SOLID PROPERTIES // 

1 60H RPS T H* GAP K SMAX , 

2 60 H THIKNESS E TO ALFAF ALFAG , 

3 10H ALFAR / 

4 60 H LB/IN N/M MINCH E-6*M L3/IN2 £9*N/M2 PSI E3*N/, 

5 60HM2 MIL E-6*M E6*PSI E10*N/M2 LB/IN N/M E-6/F /K E-6/F /, 

6 1 1HK E-6/F /< ) 

GO TO 300 

230 GOMTRP=. FALSE. 

WRITE (6,400) 

WRITE (6, DATA) 

WRITE (6, SI) 

GGMTR Y= .FALSE. 

FLUIDP=. FALSE. 

FLUID = .FALSE. 

SULIOP=. FALSE. 

SOLID = .FALSE. 

TER MUP=. TRUE. 

PRINT 231 

231 FORMAT ( 47X , 35HV AR I A T l ON OF T EMPERATU RE/ S PEED // 

1 bOH RPS T H* GAP K TFOIL , 

2 60H TGU I Dfc TROTOR VISCOSITY CONDUCTIVITY , 

3 10H DENSITY / 

4 60H LB/ IN N/M MINCH E-6*M LB/IN2 E9*N/M2 F K , 

5 60H F K F K E-8*REYN NS/M2 B/HR*FT*F J/MSK , 

6 10HLBF* S2/ I N4 ) 

300 RETURN 

401 STOP 

400 FORMAT ( IH1 ) 

END 
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SUBROUTINE EQN1 

COMMON TH1KNS, RADIUS, DG, RGUIOE ,RC,DD, SPACE, 3B,RG,R0 , 

1 T, TZ, TT, ET, TZP,TP, EPS, EPS23, V ,RP S ,Y , CAY, 

2 HSTR, H STAR , LTU, PIR2, THETA, THETA1, TH2 , 

3 CTH2, STH2 , THETAG, ThfcTGI, P i , OE G , FOA, FG A , 

4 PAIE, COMP, OENSTY, PA, EMU, CUNDCT , CF, EMOL , C,AIE, 

A CON SPR , TR VI SC » TRC QND , C V I SC , CCONO , CCND,VIS, 

3 ELTOK, RGBS , ALF , S I GBRO, S I GBR , DBDYK , DSGOB , DS3R0T ,DSGDR , 
B OB DR , 

6 THERML, ALF AF , ALFAG, ALFAR, TROOM, D TF , DTG , DTR , 

8 TROTOR , DTRCT , TROT MX, TGU IDE, TGIOMX ,TFOIL, 

9 QTHIK, THIKMX,DTZ,T ZMX , DRO , ROMX , DRG, RGMX , ORPS , R PSMX • 

7 P EXT , E, FCOEF , ENU, D ,SMAX 

LOGICAL TERMOP, FLU I DP , SOL I OP, GOMTRP ,CONSPR 

LOGICAL TERMO, FLUID, SOLID, GOMTRY ,SKIP , READRC 

COMMON TERMOP, FLU l DP, SOL IOP, GOMTRP ,ISCAN,IFOIL 

COMMON TERMO, FLUID, SOLID, GOMTRY , SKIP, READRC 

COMMON AT ( 5 ) , AX3(5), A SI GBR < 5 ) » ABSAI5I 

COMMON AFP, AGP, ARP, AFINT, AGINT, ARINT, AFINTP, 

L AGINTP , ARINTP, VtSCP, TINT, HINTP , HINT, 

2 CYINTP, CYINT, SMINTP, SMINT, ROINT, RGINT, RCINT, 

3PAINT, PAINTP, TRMINT, DNS I NT , VISINT, VSINTP, CNDINT, 

ATKIN TP, TKINT, EINT, TZINT, TFINT, TGI NT , TR I NT , EP , E I NTP , 

5 TCNV * CAYCNV, DNSCNV, VISCNV , CONCNV 
IF( V.NE.O.) GO TO 999 

T=TZ 
TP=T/ET 
TZP=TP 
999 A VG=4 . 0 

1 DO 10 1=1,25 

IF (V.NE.O. ) GO TO 5 
EPS23 =0. 

HS= 0.643 
GO TO 6 

5 EPS=V/ETU 

A I E= PAIE*EPS**2/TP 
C=COMP/TP 

EPS23=( EPS/TPJ^^O. 6666666667 

HS= 0. 643+0. 286*AIE+l.905*AIE*AIE-0. 183 /C 

6 IFID.Nt.O.) GO TO 7 
0SG3R=0. 

GO TO 8 

7 11=1 

CALL SLACK ( T , S IGBR , RADI US , RGUlDt, RC,1I) 

DSGBR= SIGBR-SIGBRO 

8 T PN= TZP + ( -2.*Y*STH2 

1 +THETA * HS*-EPS23 + D SGBR*2 . -THERML ) /(ELTOR+ PEXT) 

HSTR=R0 * HS * EPS23 
IFITPN.GT.O. ) GO TO 9 
TEMP= TP 
TP = TP/2. 

GO TO 91 
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9 IF (AHS << TPN-TP)/TPN).LT. 0.0005 ) 00 TO II 
TEMP=TP 

TP=!TPN+AVG*TP)/ (l.+AVGJ 
91 IF( IFOIL.EQ. I ) CALL TEMPF 
T- TP*ET 

10 CONTINUE 

IF ( AVG. GT. A. ) GO TO 101 
AVG=8 . 

GO TO 1 

101 PRINT 80, TEMP, TP 

80 FORMAT! 35H NO CONVERGENCE AFTER 50 STEPS, 6HTP0L0= , 

l E16.8, 6H TPNE W= E16.8) 

11 CONTINUE 

HSTAR ■= HSTR*1.E 6 
IF! D.NE.O.) GO TO 27 
0S8R DT=0 . 

0 SGQB=0. 

DSGDR=0. 

GO TO 28 

27 XXX-= ( ( AS I GBR ! 1 ) - ASIGBR (3 ))/ (AT l.l) -AT!3)> 

1 ( ASIG8R(2) — ASIGBR! 3) ) / ! AT C 2)— A T ( 3 ) ) J / !AT<1) -AT!21) 

YYV = ( ( ASI GBR ! 1 ) -ASIGBR!3))/ !AT!1)**2 - AT(3)**2) - 

1 ! AS IGtiR 12 )- AS I GBR l 3 ) ) / (AT(2)**2 - AT!3>**2)) 

2 * ! AT! 1 ) + AT ( 3 ) J * < AT ! 2 > + AT! 3) ) / (AT!2) - AT l 1 ) ) 

0S8RUT =2. * XXX * T + YYY 

T=TP*ET 
II =4 

CALL SLACK l T, SIGBR, RADIUS, RGUIDE, RC.II) 

28 CAY = 3.*STH2*((2. *STH2 - DSGDB * DBDYKI / RADIUS ) / 

1 l (ELTOR+PEXT)/ET - EPS23 * ( {-. 286*A I E -3.81*AIE**2 - 0.183/C 

2 -0.66ob 67*HS) / T ) * ( THETA+2 . *( DSGOR+ OSGDB*QBDR ) )-2 .* DSBKDT ) 
SM AX = T/THIKNS*- ET/ ! 2 . *RGU I Dt ) 

SMI NTP= SMAX*6. 3947E-5 

C SMINTP=SMAX*6. 3947E3 *(l.E-3 FOR PRINTING) 

TI NT=TCNV*T 
HIi\ITP=HSTAR*2. 54E-2 
HINIT=HINTP*l.E-6 
CYINTP=CAY*CAYCNV 
CYINT = CYINTP*1.E9 
R0UN0=100 . 

IFtCAY. GT. 10000. J GO TO 13 
ROUNO= 10. 

13 CAY- ROUND* I F I X ! ( C A Y+RUUNO/ 2 . ) /ROUND ) 

C IF1SKIP) LEAVE A BLANK OUTPUT LINE 

IF! .NOT. SKIP) GO TO 15 
SKI P=. FALSE. 

WRITE (6,14) 

14 FORMAT ( 1H ) 

15 CONTINUE 

IF(GOMTRY) GO TO 200 

IF(FLUID) GO TO 210 

IF! SOLID ) GO TO 220 
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n l TtRMO) GG TO 230 

200 PRINT 201 » RPS, T, TINT , HS, HSTAR , HINTP, CAY, CYINTP, SMAX, 

I SMINTP, ROtROINT, RG, RGINT , RC, RCINT , THETA1, THETG1 

201 FORMAT! IX, F6.0, F6.2,F6.0 , F6. 3 , F6. 0 , F5 . I , F7 . 0 , F6 . 3 , 1 X , 

l F8.0,F5.1, F8.4, F7.4, F8.4, F7 . 4, F 8 .4, F 7.4 , F 8. 2 , F 7. 2 ) 

GO TO 300 

210 CON 01= CONDCT/CF 

PRINT 2 L 1 ,RPS , T , TINT, HS, HSTAR, HINTP, CAY, CYINTP, 

1 PA, PAINTP, TROOM, TRMINT, EMOL, OENSTY, DNS I NT , VISCP, 

2 VSINTP, CONDI, CNOINT 

2 1 1 FORMAT ( IX,F6.0,F6.2,F6.0,F6.3,F6.0,F5.1, F7.0,F6.3, IX , 

IF6.I, F6.I, F5.L, F6.1, F5.1, EI0.3, F7.2, F7.3, F7.3, EI0.3, 
2 E10.3) 

GO TO 300 

220 PRINT 221, RPS, T, TINT, HS, HSTAR, HINTP, CAY, CYINTP, 

1 SNA X , SMINTP, TT , TKINT P, EP, EINTP, TZ,TZINT, 

2 AFP , AFINTP, AGP , AGINTP, ARP, ARINTP 

221 FORMAT! I X , F6 .0 , F6. 2 , F6.0, F6 . 3, F6.0, FS.l, F7.0, F6.3, IX , 

1 F 3 . 0 , F 5. 1 , F5.2,F5.1, F7.1, F7.2, F3.1, F6.1, F5.1, F5.I, 
2Fb . I , F5.I, F6.1, F5.1 ) 

GO TO 300 

230 COND 1= CGNOCT/CF 

PRINT 231, RPS, T, T INT, HS, HSTAR, HINTP, CAY, CYINTP, 

1 TFOIL, TF I NT , TGUIOE, TGINT, TROTOR, TRINT, 

2 VISCP, VSINTP, C0N01, CNOINT, DENSTY 

231 FORMAT ( IX ,F6.0, F6.2, F6.0, F6.3, F6.0,F5.1, F7.0, F6.3, IX, 

1 F 6. 1 , F 6. 1 , F6.I, F6.1, F6.1, F6.I, F7.3, F7.3, E10.3.E10.3, 

2 E10.3) 

300 RETURN 

END 
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.SUBROUTINE SLACK ( TX, S1G6RX.RQX, RGX,RCX,I) 

“COMMON THIKNS.RAOIUS. OG.RGU IDE , RC, OD, S P ACE , BB, R G, RO , 

1 T , T Z, TT, ET, TZP,TP, EPS , EPS23 , V , RP S , Y , CAY, 

2 HST R , HSTAR, ETU, PIR2, THETA, THETA1, TH2 , 

3 CTH2, STH2 , THETAG, THETG1, P I , DEG, FOA, FGA, 

4 P A I E , COMP, DENSTY, PA, EMU, CONDCT ,CF,EMOL,CZ,AIE, 

A CONSPR,TRVISC,TRCONO,CV ISC.CCONO, COND,VIS, 

5 ELTOR, RGBB , ALF , S I GBRO , S I GBR , DBDYK, DSGOB, OSBROT , DSGDR , 
B OBOR, 

6 THERML , ALFAF, ALFAG, ALFAR, TROOM, DTF, DIG, OTR, 

8 TROTOR »OTRGT ,TROTMX»TGUIOE, TGIDMX ,TFOIL, 

9 DTHIK, THIKMX,DTZ,TZMX,ORO,ROMX,DRG,RGMX,DRPS,RPSMX, 

7 PE X T, E, FCOEF, 6NU, 0 ,SMAX 

LOGICAL TERMOP, F LU I OP , SQL I DP , GOMTRP ,CONSPR 

LOGICAL TERMO, FLUID, SOLID, GOMTRY ,SKIP ,READRC 

COMMON TERMOP, FLU I DP , SOL I OP , GOMTRP , ISCAN, I FOIL 

COMMON TERMU, FLUID, SOLID, GOMTRY , SK i P , RE ADRC 

COMMON AT ( 5 > , AX3I5), AS I GBR 15), ABSAI5) 

COMMON AFP, AGP, ARP, AFINT, AGINT, ARINT, AFINTP, 

1 AGINT?, ARINTP, V I SCP , TINT, HINTP, HINT, 

2 CYINTP, CYINT, SMINTP, SMINT, ROINT, RGINT, RCINT, 

3PAINT, PAINTP, TRMINT, ONSINT, VISINT, VSINTP, CNDINT, 

4TKINTP, TKINT, E I NT , TZINT, TFINT , TGI NT , TR I NT , EP, E I NTP , 

5 TCNV, CAYCNV, DNSCNV, viscnv ,CQNCNV 

DOUBLE PRECISION A, A2 , C, SH3 , CH3, TH3, PARTI, PAR T2 , P,Q, 

1 X.I ,X1A,X2,X2A, X3 , 8SA,EPSA,TA,SHI,CH1,SH2,CH2, 

2 TERM2, TERM4, TERM5 , TERM7 , S I 8RE3 ,S I GM 

EQUIVALENCE (X1,X1A), (X2.X2A) ,( SI BR E3 , S IGM ) 

C WRITE (6,9) 


9 FORMAT ( IH , 

3X , 7HT EN S 1 ON , 5 X , 3HEP S , 

6X, 

2HF0, 

6X , 

2HFG, 

4X , 

1 4HSCI/3, 4X , 

2 2H J , 

4HSG/B, IX , 7 HB/ ( R *E ) 

IX, 7HTENSI0N,5X,3HEPS , 

6X , 

2 HFO , 

6X , 

2 HFG, 

4X , 


3 4HS0/3, 4X, 4HSG/B, IX, 7H8/ ( R*E ) ) 

SBX = SQRT ( RCX**2 - (ROX + RGX ) ** 2 ) 

I F ( I . N E . 4) GO TO I 
ROBB = 1.10* R0X/86X 
R08RG = 1.10* ROX/RGX 
BBYRQ = BBX / <ROX*l.IO) 

R0BRG2 = R03RG**2 
RGBRO = 1./ ROBRG 
X3 = AX 3(3) 

GO TO 3 

1 ROBB = ROX / 3 B X 
R08RG= ROX/RGX 
BB YRO= 8BX/R0X 
ROBRG2 - R03RG**2 
RGRRO = RGX/ROX 
TA = TX +0.35 

C ESTIMATED FIRST GUESS 

2 X3 = ( RGX* FGA+83X+R0X* FO A ) *DS QRT < T A/ D > 

3 J = 0 

IF(X3.LT .85. ) GU TU 110 
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SH3=1 .E43 
CH3-1.E43 
T H3= 1 . 

A=-l. 

A 2=1. 

C=RO iiKG* ( X3— 1 . ) 

GO TO 120 
110 CH3=0CGSH< X3) 

SH3 = DS I NH ( X 3 ) 

TH 3 =0TANH( X3 ) 

A= - ROBRG/ SH3-1./TH3 
A2 = A **2 

C= ROBRG* X3-R0BRG/TH3 -1./SH3 
120 P= R0BRG*X3-(R0BRG-1. )*OTANH( X3/2- > 

Q = 2.+ R0BRG*(X3**2+2. ) -2.*X3 *(1./SH3+ ROBRG/ DTANH ( X3 ) 
PARTI = P*ROBRG 

PART2= Q*( R0BRG2+R0BRG) -P**2*R08RG 

IF { PART2.GE.0) GO TO 5 
PART 2 = 0. 

J=1 

5 PART 2 =OSQRT ( PART2 ) 

X2 A= (PARTI + P AR T2 ) / ( R0BRG2 + ROBRG) 

X1A= P-R0BRG*X2A 

BSA= X2A - X 1 A 

EPS A= B6YR0/BSA 

FGA= R0BRG*EPS A* ( X3 -X2A) 

FO A= X 1 A* EPSA 
F0=F0A*DEG 
FG = FGA*DEG 

C SOOBA = SO/B SOLUTION A ETC 

SOBBA = X 1 A* EPS A *R0BB 
SGBBA= ( X3 -X2A) *EPSA*ROBB 
BBRO E A= B3YR0/EPSA 
SHI =DS I NH ( X l ) 

CHI =OCOSH(X1) 

SH2 =QS I NH ( X2 ) 

CH2 =DC0SH(X2) 

I F ( X3 .LT .8 5 . ) GO TO 210 
219 TERM2 = - R08RG*DEXP(X2— X3)«- DEXP(-X2) 

TERM 4= 2 . *R OBRG * ( -R0BRG*DE XP( X2-X3 ) -DE XP (-X2 ) ) 

IF (XI .LT.O. ) GO TO 220 

TERM 7=— 2 . * ( ROBRG*DE XP( X l— X3 ) + OEXP(-Xl)) 

1 — ( 0 . 666667* XI + A)*X1**2 
GO TO 225 

220 TERM7=-Xl**3/3. *2.*<A+X1> 

GO TO 225 

210 TERM2 = (— R0GRG*SH2 + DS I NH ( X3-X2 ) ) / S H3 

TERM4= 2.*R03RG* {-R0BRG*CH2 - OCOSH( X3-X2 ) ) /SH3 
IFtXl.LT.O.) GO TO 240 

TERM7=-2.*(R0BRG*CH1 +DC0SH(X3-X1 ) )/SH3- 
1( 0 . 66666 7 *X 1 +At*Xl**2 

GO TO 225 
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240 TERM7= -Xl**3/3.*2.*(A+Xl) 

225 IF(X3.GT.16.)G0 TO 22i 

TERM5 = (A2+l.l*DSINH(2.*X3)/4.+A*DC0SH( 2.*X3)/2. 

GO TO 228 

221 TERM5 = 0.5*R0BRG2 

228 SIGM =(-A/2.+< A+Xl )**2*X2+ ROBRG2’ !t (X3**3-X2**3)/3. 

1 +( A2-1. )*X3/2.+C**Z*l X3— X2)-R0ijRG*C*<X3**2-X2**2) 

2 +T ERM4+TERM5 +TERM7J/2. 

SIR2B3 = SIBRE3 /BSA**3 

S IG8RX= SIR283/R0BB**3 

TA= ( 6SA/8BX)**2*D 

AT( 11= TA 

AX3( I ) = X3 

ASI GBR ( I )= S1GBRX 

ABSA ( I ) = B SA 

GO TO (2 4,25,26,34,35 I 

24 1=2 

TA= TA- 0. 70 
GO TO 2 

25 X3 = A X3 ( L ) *■ (AX3 12)- AX3 ( 1 ) )*( TX-AT < 1 ) ) / C AT C 2 ) —AT C 1 » ) 

1 = 3 

GO TO 3 

34 X3 = AX3(3>*1.10 

DSGOB= ( AS I GBR ( 2 )— AS I GBR ( 1 > ) / { ( ABSA { 2 )- A6 SA ( 1 ) ) *EP SA ) 

1=5 

GO TO 3 

35 BS I GBR = ASIGBK(4) + (ASIGBR(5)— AS I GBk 14))* (ABSAI3) -ABSA (4)) / 

1 ( ABSA ( 5) — AB SA ( 4 ) ) 

nSGOR=(0SIGBR*l.lO - ASI GBR ( 3 ) ) / 0.1 

26 RETURN 
cND 
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SUBROUTINE TEMPF 

C DM MON THIKNS, RADIUS, OG,RGUIOE,RC,OD, SPACE, BB,RG,RO , 

1 T, TZ, TT, ET, TZP , TP , EPS,EPS23,V,RPS,Y, CAY, 

2 HSTR, HSTAR, ETU, PIR2, THETA, THETA1, TH2, 

3 CTH2, STH2 , THETAG, THETGI, PI ,DEG, FOA, EGA, 

4 P A I E, COMP, OENSTY, PA, EMU, CUNDCT , C F , EMOL , C,AIE, 

A CONSPR.TRVI SC, TRCONO,CV ISC, CCOND, CONO.VIS, 

5 ECTOR, RGBB , ALF , S l GBRO , S I GBR , DBDYK , DSGDB, DSBRDT , DSGDR , 

B OBOR, 

6 THERML, ALFAF, ALFAG, ALFAR, TROOM, OTF, DTG, DTR, 

8 TROTOR » DTROT » TROTMX, TGU IDE » TGIOMX ,TFOIL, 

9 DTrtlK, T H I KMX » OT Z » T ZMX , ORO »ROMX » DRG, RGM X, DRPS, RP SMX * 

7 PEXT, E, FCOEF, EiMU , D , SMAX 

LOGICAL TERMOP, FLU I OP , SOL I DP , GOMTRP ,CONSPR 

LOGICAL TtRMO, FLUID, SOLID, GQMTRY .SKIP , R E A DR C 

COMMON TERMOP, FLUIDP, SOLIDP, GOMTRP ,ISCAN,IFOIL 

COMMON TERMO, FLUID, SOLID, GOMTRY ,SKIP,R£ADRC 

COMMON AT I 5) , AX3I5), AS I GBR ( 5 ) » ABSAI5) 

COMMON AFP, AGP, ARP, AFINT, AGINT, ARINT, AFINTP, 

1 AGINTP, ARINTP, VISCP. TINT, HINT P » HINT, 

2 CYINTP, CYINT, SMINTP, SMINT, RG I NT , RGINT, RCINT, 

3 PA I NT , PAINTP » TRMINT, DNS I NT , VISINT, VSINTP, CNDINT, 

4TKINTP, TKINT, EINT, TZINT, TFINT, TGINT, TR INT , EP , E INTP , 

5 TCN V , CAYCNV, ONSCNV, VISCNV .CONCNV 

I F ( CONS PR ) GO TO L21 

EMU = VI S* I TRVI SC+460.+CVI SC ) * ( ( TROOM+ ( DTR+DTF )/2 .+460. ) / ( 
2TRVISC+460. ))**li5/ ( TROOM +IDTF +UTR)/2. +460. + CVISC) 

CONUCT= CONO* CF* ( TRC0N0+460 . +CCOND )* II TROOM + ( DTF + DTR ) /2 . +460.)/ 
1 I TRC0N0+460. ) ) **l . 5/ { TROOM+ ( OTF+DTR ) /2 . +460. +CCOND ) 

GO TO 122 

121 EMU= VI S 
CGNDCT=CQND*CF 

122 CONTINUE 

IF( IFOIL.NE. 1) GO TO 10 

T FO I L - ( |£MU*V**2/ (2.*C0N0CT > +TROTOR) * I DD-HSTR) /DD 

L + TGUIDE* HSTR / DD 
10 DTF = TFO I L-TROOM 
C OENSTY = LBF*SEC**2/ IN**4 

ETU=ET/(6.0*EMU) 

DENSTY = IPA/I 1544. *(460.+ TROOM+ (DTF + DTR ) /2 . ) / E MOL ) ) 

1 / I 32 . 2 *144.) 

P AI E = DENSTY *E T *RADIUS/(72. *EMU**2) 

5 THERML = ( AL FA F* AL F *DT F 

1 + 2.*ALFAF * (RGUIDE *THET AG) *DT G 

2 -ALFAG*2.*38*DTG -ALFAR* THETA* RAD I US *DTR -2. *AL FAG* THETAG 

4 *RGU I DE*DT G) /RADIUS 
TF I NT=( TFUI L+459.67) /I .8 
DNS I NT= DENSTY* ONSCNV 
CNDI NT=CONDCT*CONCNV 

VI SC P = E MU* 1 . E 8 

VS INT P = EMU*V I SCNV 

VI ST NT= VSINTP* L.E-4 

RETURN 

END 
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SYMBOLS 


distance between origin "O" and line connecting the centers 
of the guides. Fig. 3 (k refers to kth foil sector, subscript 
k is omitted for axisymmetric case), in. (m) 

2 2 

function defined in Eq. (20b) lb/in (N/m ); also, integration 
constant Eq. (A19) 

points defined in Fig. A1 

distance between centers of guides, Fig. 3 (k refers to kth 
foil sector; subscript k is omitted in axisymmetric case) , 
in. (m) 

partial foil length defined in Fig. 3 and Fig. A1 (k refers to 
the kth foil sector, i and e refer to inlet and exit respectively), 
in. (m) 

integration constant, Eq. (A20) 

Sutherland constant for thermal conductivity, Eq. (30), 

°F (°K) 

Sutherland constant for viscosity, Eq. (29), °F (°K) 
foil thickness, in. (m) 

distance defined in Fig. 2 for reference state, in. (m) 

flexural rigidity of foil = Ed 3 /|2 (1- V 2 )]/ lb/in (N/m) 

2 

Young's modulus, psi (N/m ); also integration constant Eq. (A21) 
subscript referring to the exit region 

components of resultant force of foil support system on the 
rotor, lb/in (N/m) 

subscript referring to guide or shield- damper, Fig. 2 

distance of rotor to kth shield-damper. Fig. 2 (subscript k 
is omitted in axisymmetric case), in. (m) 

distance of rotor to shield-damper in reference state, in. (m) 



h , h clearance in uniformity region of kth foil sector (subscript 

k is omitted for axisymmetric case), in. (m) 

H*,H* dimensionless clearance in uniformity region of kth foil 

sector (subscript k omitted for axisymmetric case) 


i 

k 


«e: 

M 

O 



r 

c 


r 

g 

r 

o 

R 


s 

S 1' S 2' S 3 

T 

o 



w 


W 

U 


subscript referring to the inlet region 

bearing stiffness per unit width lb/in^ (N/m^); also subscript 
referring to kth foil sector 

total foil length between clamping points on foil guides, 
referred to kth foil sector, in. (m) 

magnitude of foil length 1 in the reference state, in. (m) 

)>► 

effective foil length defined by Eq. (11), in. (m) 

molecular weight of lubricating fluid 

position of rotor center in reference state. Fig. 2 

pressure of fluid entering the bearing, absolute, psi (N/m ) 

local radius of curvature of foil , in. (m) 

local radius of curvature of perfectly flexible foil, in. (m) 

distance of center of guide from origin O, Fig. 3, in. (m) 

radius of foil guide, Fig. 3, in. (m) 

radius of rotor in reference state, in. (m) 

universal gas constant, in lb/lb-mole °F (mN/kg-mole °K) 

length along foil, in. (m) 

values of s at points k^,k^,k^, Fig. A1 , in. (m) 

initial tension per unit width of foil (preload), lb/in (N/m) 

tension per unit width of kth foil segment (subscript k may 
be emitted for axisymmetric case), lb/in (N/m) 

deviation of finite stiffness foil contour from flexible foil 
contour, in. (m) 

w/(r Q £ 2 ) 

surface velocity of rotor, in/sec (m/s) 
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x,y 


V y k 


^f/fk 


components of displacement of rotor center. Fig. 2, in (m) 

component of displacement of rotor center in the kth auxiliary 
coordinate system, Fig. 2, in. (m) 

coordinate normal to foil, Fig. 4, in. (m) 

coefficient of thermal expansion of kth foil, (subscript k 
is omitted in axisymmetric case), l/°F (l/°K) 


oC 


°^ik' °^ek 


coefficient of thermal expansion of kth shield-damper-guide, 
(subscript k is omitted in axisymmetric case), l/°F (l/°K) 

coefficient of thermal expansion of rotor, l/°F (l/°K) 

inlet and exit angles defined in Fig. 3, referred to kth 
foil sector, rad (rad) 


Pik' hk 

Kk 

«k 


inlet and exit angles defined in Fig. 3, referred to the 
kth foil sector, rad (rad) 

change in a relative to reference state, in. (m) 

G K. 

elongation of kth foil “ £ - £ , in.(m) 

K O 

angular position of kth foil sector referred to the foil 
sector (k = 1) , Fig. 2, rad (rad) 


£ yW >- 0 2 ) 

<£b, £b wP-’ik '<S b . change in b ,b 

C " state, in. (m) cK 




with respect to reference 



equivalent length change of kth foil sector defined by 
Eq. (10), in. (m) 





<Tt 

^fk 


perturbation in rotor radius (thermal and centrifugal) , 
in. (m) 

perturbation in guide radius, referred to the kth foil sector, 
(subscript k is omitted for axisymmetric case), in. (m) 

difference between value of 6”ik/Gek or < *< anc * th® va l ue 
<j 0 (reference state), in (m) 

excess temperature over reference, X - X , °F (°K) 

o 

- X (subscript k is omitted for axisymmetric case), 
tk o 

°F (°K) 
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© 

g 



(subscript k is omitted for axisymmetric case) , 
°F (°K) 



(subscript k is omittea ior axisymmetric case) , 
°F (°K) 


angles defined in Fig. Al, rad (rad) 


nominal angle of wrap in reference state, rad (rad) 

nominal angle of wrap around guide in reference state, rad (rad) 


thermal conductivity of lubricating fluid at temperature T , 
BTU/ft hr °F) (j/(m.s °K) ) a 


k, thermal conductivity of lubricating fluid at base tempera- 

" ture T b , BTU/tt.hr °F) (j/(m.s ok) ) 

p viscosity of lubricating fluid at temperature T. , lb sec/in^ 

(N.S/m^) 

n 

p., viscosity of lubricating fluid at base temperature ^ , lb sec/in 

(N . S/m 2 ) 

T? Poisson's ratio of foil material 

f density of lubricating fluid at pressure p and tempera- 

c ture X , lb*secVin^ (kg/m^) J 

cl 

< r excess of length of finite stiffness foil over a perfectly 

flexible foil defined in Eq. (A25), in. (m) 

cr ik ,cT ek value of<J* referred to inlet and exit branch of kth foil sector 

respectively, in. (m) 

X , X temperature of working fluid; assumed; X = (X + X f , )/2, (sub- 

3 Q.K 3. r IK 

script k is omitted for axisymmetric case. The subscript 
(abs) is added when absolute temperature is used). 

°F (°K) 

^ s/(r D £ ) 

3l*| 2^3 values of | at points A^ ,A^ ,A^ , Fig. A1 
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